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Abstract 

The effect of retardation on dynamical mass generation in QED in two space 
dimensions at finite temperature is studied, in the imaginary time formalism. The 
photon polarization tensor is evaluated to leading order in (where A'^ is the 
number of flavours), and simple approximate closed form expressions are found for 
the fully retarded longitudinal and transverse propagators, which have the correct 
T — > limit. The resulting Schwinger-Dyson equation for the fermion mass (at 
order has an infrared divergence associated with the contribution of the 

tranverse photon propagator; only the longitudinal contribution is retained, as in 
earlier treatments. For solutions in which the mass is a constant, it is found that 
retardation reduces the value of the parameter r (the ratio of twice the mass to 
the critical temperature) from about 10 to about 6, in agreement with a similar 
calculation in the real time formalism. The gap equation is then solved allowing 
the mass to depend on frequency (but not momentum), thus extending the study 
of retardation to the variable mass case for the first time. Solutions for T 7^ are 
obtained which join on smoothly to the correct T = solution. It is found that 
there is a critical number of flavours, A^ci^bove which no mass is generated. The 
phase boundary in the N — T plane is calculated, and agrees qualitatively with 
that found in other variable mass (but non-retarded) calculations. The r value 
remains close to 6. Possibilities for including the transverse photon propagator are 
discussed. 
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1. Introduction 



The study of quantum electrodynamics in (2+1) dimensions is of con- 
siderable interest, due to its possible relevance to long- wave length models 
of 2-D condensed matter systems, particularly those which might apply to 
high Tc superconductors [13], [14]. QED^ may also give us insight into the 
phenomenon of dynamical symmety breaking in theories such as QCD. At 
zero temperature, a considerable amount of work has been done [1] — [8]. 
In Appelquist et al [2] it was shown that, using a 1/N expansion in the 
Schwinger-Dyson (S-D) equations, there existed a value Nc — 32/7r^ above 
which no fermion mass was dynamically generated, N denoting the number 
of fermion flavours in the theory. This was found to be qualitatively still true 
to order 1/iV^, changing by a factor of 4/3. [5]. In Pennington and collab- 
orators [3], by contrast, no N,, was found; instead, the dynamically generated 
mass fell exponentially with increasing N . This work adopted a more general 
non-perturbative approach to the S-D equations. On the other hand, an al- 
ternative non-perturbative study by Atkinson et al [6] claimed that there was 
indeed a critical number of flavours. More recently, Maris [7] has confirmed 
the existence of an with a value of about 3.3,more or less independent 
of the choice of vertex ansatz, by considering the coupled S-D equations for 
the photon and fermion propagators; and Kondo [8] has studied the problem 
using the non-local gauge technique. 

It is important to extend this type of analysis to finite temperature, because 
of its possible relevance to physical applications. Here the significant param- 
eters are Tc and r, where r is the ratio of twice the zero temperature fermion 
mass to the critical temperature, Tc-T^ is the temperature above which chiral 
symmetry is restored, and at which there is a phase transition from the su- 
perconducting phase to the normal phase in the model discussed in [13]. In 
Dorey and Mavromatos [9] a calculation was done using the S-D equations, 
with the Matsubara finite-temperature formalism, in which the fermion mass 
S was taken to be constant, and only the Aqo component of the photon prop- 
agator was used in the instantaneous approximation-that is, all frequency 
dependence in the photon propagator was neglected. In this calculation the 
authors found a value of r ~ 10. If we identify the fermion mass with the 
order parameter in B.C.S- like theories, this value is much larger than a typ- 
ical B.CS value, which is r ~ 3.5. In Aitchison et al [10] the constant mass 
approximation was relaxed, and instead a momentum -dependent solution 
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E(T', p) was calculated. However once again only Aqo in the instantaneous 
approximation was retained, and it was found that r ~ 10. The r value was 
therefore insensitive to this refinement in the calulation, although the values 
of Tcks and S(T = 0, p = 0) had changed considerably. It was shown in 
Aitchison and Klein [11] that this value of r also survives the inclusion of a 
form of wavefunction renormalization. In both the calculations of [10] and 
[11] the critical value of N was ^ 2. It is important to stress that no such 
Nc existed in the constant E case [9]. 

One obvious problem with the instantaneous approximation in the Mat- 
subara formalism, is that it cannot reproduce the well-studied zero-temperature 
limit. This is because, as T — 0, all frequency components should be in- 
cluded, while the instaneous approximation retains only the n = one. To 
avoid dealing with a large number of frequency components, the first calcula- 
tion to include retardation in the S-D equations for QED^ at finite tempera- 
ture [12] used a real-time formahsm, in which the correct T — > was ensured. 
Making the constant mass approximation, it was found that r was signifi- 
cantly reduced to a value of about 6, from the non-retarded value of about 10. 

The calculations of [12], however, did not retain the exact expressions for 
the longitudinal and transverse photon self-energies,because of their awkward 
behaviour near zero three-momentum ( the amplitudes are non-analytic at 
the origin, at finite temperature). Instead, a variety of simpler "average" 
self-energies were used in order to simplify the calculation and to be able 
to compare the results more easily with those of [9]. In the imaginary time 
formalism, this difficulty concerning the non-analytic behaviour near zero mo- 
mentum does not arise, and the exact self-enegies can be employed. Our first 
aim in this paper, therefore, is to include retardation effects, in the constant 
mass approximation, using the imaginary time formalism. In Section 2 we 
restate for convenience the results given in [13] for II^i,, the photon polariza- 
tion tensor to leading order in 1/N. These results involve certain integrals 
which we first evaluate numerically, and then find simple approximate closed 
form expressions which retain the correct T — ^ limit. In section 3 we shall 
then extend the (constant mass) work done in [9] by using a fully retarded 
longitudinal propagator. We compare the results with [12], and find rather 
close agreement. 

In section 4 we shall go further and allow the fermion mass to depend 



2 



on the (discrete) frquency,still retaining the fully retarded longitudinal self- 
energy. We believe this is the first time that retardation has been introduced 
into a "non-constant-mass" S-D calculation at finite temperature. Although 
we cannot compare our results directly with [10], there are interesting fea- 
tures which remain the same, namely A^c ~ 2 and apo/a behaviour similar to 
that of the \p\/a behaviour of the momentum dependent solution. Our val- 
ues of S(po = 0) are found to be in rough agreement with those of S(p = 0) 
in [10]. It would, of course,be interesting to include a dependence of S on 
momentum as well as on frequency, but that is a much harder problem in 
the retarded than in the instantaneous case. 

In section 5 we shall discuss ways of introducing a transverse contribution 
into our mass-gap equation, which we neglected due to a logarithmic i-r 
divergence in the zeroth mode. We shall also discuss other possible extensions 
to our calculation, and we conclude by looking at the plausiblity of our 
calculations in the context of QED^ as a model of superconductivity. 

2. An approximate form for the full photon propagator to leading 

order in 1/A^ 

The Lagrangian of massless QEDs with N flavours is 

L = -i/Af,,r +Y.U^^ - 4yi'^ (1) 

i 

where a^^ is the vector potential and i = 1, 2... A". We have also chosen a re- 
ducible representation of four-dimensional matrices for the Dirac algebra. 
Due to the choice of representation eqn (1) has continuous chiral symmetry 
as discussed in [2]. 

Following [13] we calculate the full photon propagator in this theory to 
leading order in 1/A" . We then look for simple closed form expressions to 
approximate those integrals in the calculation which cannot be evaluated 
analytically. In these calculations we shall be working in the Landau gauge, 
in Euclidean space. We first note that the most general photon propagator 
in the Landau gauge must take the form: 




where is the longitudinal projection operator and B^i, is the transverse 
projection operator. To leading order in 1/N,IIa and Ub are the contribu- 
tions from the loop diagrams shown in fig. 1 and are related to the polarization 
tensor, U^i, ,at order The projection operators take the following forms: 

and = 5f,i (^ij - ^^^j (3) 

By using the properties of the projection operators it is easy to show that 
the inverse propagator must be: 

= (p' + n^(p)) + (p' + \Ib{p))b^,{p). (4) 

^From this expression it is easy to relate 11^^ to 11^ and by summing the 
diagrams in fig. 1 : 

n^, = Ha V + ^bB^,. (5) 

Then using the explicit forms of A^^i, and B^i, in (3) with the additional 
requirement that PfjU.^" — one obtains the following forms of 11^ and Ub- 

Ha = noo4 

Ub = Uii-Uoo^. (6) 

We now proceed with the calculation of H^j, at finite temperature using the 
Matsubara formalism. We work in imaginary time in which both fermonic 
and bosonic frequencies are discrete. Fermionic frequencies have the form 
Pof = "^ijf^ + 1/2) and bosonic frequencies have the form po6 = where 
m is an integer. From now on we shall denote the modulus of bosonic 3- 
momenta by pb and that of fermionic 3-momenta by p/.In our calculation we 
need only consider the elements Hoo and Iljj. Hoo and Iljj are calculated to 
be [13]: 

Hoo = n3-4ni- 

Pb 




n,,- = ii,-\5i^-^\+iiAj (7) 
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where 



and 



2% Jo ^ Drn{x,Ph,p) 

2/3 Jo Dm[x,pb,l3) 
a 

U3{ph,P,m) = — / dxln{4:Dm{x,pb,(3)) 
np Jo 



Djn{x,Pb, P) — cosh^{Ppb^x{l — x)/2) — sm^(xm7r). 



(8) 



(9) 



m denotes explict frequency dependence in the above expressions. 
Using the result 



Ha = Hi + 

Pob 



(10) 



we find that = na,!!^ = 111 + 112. Alternatively one can obtain the above 
expressions for Ha and Hb by neglecting 112 and setting Hi = lis, which is 
found to be a good approximation for Iljj and IIoo when m 7^ 0. For m = 
one can find these expressions without again using eqn (10), by noticing that 
n2(m = 0) = and po6 = 0. 

In our treatment of Hi, 112 and Ha it will be convenient to consider the 

m = mode separately. For m = we need only consider ITi and Ha, since 
n2(m = 0) = 0. For Tl^{m = 0), which we shall denote by Ilg, an accurate 
closed form expression has been given in [10], namely 



n° = 1/8 



a 



, 1^ 161n2 

IpIp h — z — 



TT 



TT 
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(11) 



We are able to find a similar expression for 11° — H^ijn — 0) by noting that 



n? = (|p|) 



which gives, on combining eqn (11) and eqn (12), 



«/^|p|/5' 



exp 



71 
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(12) 



(13) 



Wc now evaluate the integrals Ili{pb, P, m), Y[2{j>h, P, ™) and Hsi^Pb, P, fn) for 
m 7^ 0. We show the numerical results of these as functions of pi,P and m in 
figs 2a, 2b and 2c. We see that both IIi and lis vary little with the explicit 
frequency index m. Also we are able to show that < 112 < (^^) (^)- 
and Ha are well approximated by 

Hi = Ha ^ 1/8 (^^ j (14) 

for m 7^ as figs 3a and 3b show. Since 112 is bounded by (^^^ it is 

always smaller than IIi for m ^ 0, since IIi > (^'^^ (^^^ ,so we shall neglect 
it in n^. 

We now have two expressions for the full photon propagator: 

A., -^A^ + ^l^^iorpoby^O. (15) 
Pb + c^Pb/^ P^ + ocpbl^ 

The structure of our propagator is consistent with the zero temperature re- 
sult, to order 1/N: Ha — ^b — ap/S. One should note that although in the 

zero temperature limit the theory is Lorenz invariant, Lorenz invariance is 
clearly broken at T 7^ 0, as (15) shows. This is, of course, because a preferred 
frame of reference is provided by the heat bath. 



3 The Schwinger-Dyson equation and its solution for constant E 

The full Schwinger-Dyson equation for the fermion propagator at non- 
zero temperature ksT — 1/P is given by 

s?'iPf)-sp-'{pf)-^ Y: ^,rsF{kf)A,^{kf-pf)r^{kf-pf,kf). 

(16) 

We now truncate eqn (16) by working to leading order in 1/N , in which case 
r*^ is replaced by its bare value 67"^ and we use our form for A^j,, the full 
photon propagator, given in (15). As before [9] we neglect wave function 
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renormalization; wc comment further on this in the next section. On taking 
the trace of eqn (16) we get the following equation for the mass-gap function 
E 

a ~ /■ c^2^ E(A;/) 
S(P/) = ^__Ej(^A„(*=,-p,)j^^i^ (17) 

where 

A,,(g. = %-P/) = |q|2 + n0(|q|) + |qp + n0(|q|)'f°^^°'' = 

and a — Ne^. Now we make the approximation that E is frequency as well 
as momentum independent. We can then remove a factor of E from each 
side of the equation, and write the argument of A^^ in (17) as (^, k). We 
now do the angular integration, which is trivial. Rearranging terms gives 

[SL{a,s) + ST{a,s)]. (19) 



2iV7r 

St is the transverse contribution which is expressed as: 
roo ^ 1 1 



roo 

Sria, s) = X dx 
Jo 



x^ + pmi X^ + TT^ + a?s^ 



^ ^1 (x2 + (27rm)2 + 0.125a(x2 + (27rm)2)V2) 

( 1 1 \ 



^x2 + (27r(m + l/2))2 + a2s2 + (27r(m - 1/2))^ + a^s^ 



where x — /3|k|,a = q;/3, and s = T,/a. Sl is the longitudinal contribu- 
tion and is exactly the same, except for Il'j' being replaced by Ilg. It is now 
important to notice that the integral for the 0th transverse mode is loga- 
rithmically i-r divergent. This can be seen by observing that as /3|k| — > 0, 

^1 ~^ (f) ( i28in2 ) I'^l^/^^- Proin now On, WC shall retain only the longitudi- 
nal mode; in our Conclusion we shall discuss ways of including the transverse 
mode in our calculations. 
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So we consider the following equation 

l = ^S.«...). (21) 

The integral for the m ^ modes can be done analytically. We introduce 
functions of the form: 

l{d,a,C) = J^ (3.2 + ^2) +^(^2 + ^2)1/2 3.2 + ^2 ^ ' 

Doing the integral on the R.H.S of eqn (22) gives us the following closed form 
expressions for I{d, a, c). 

= 2(a2 + Lci2) ^"((dT^) 

a /(c2__rf2)l^\ 

+ (c2-d2)l/2(a2 + c2-ci2)^'''*^''V ^ j' 

With the use of these functions we are able to write ^'^(a, s) as 

SL{a,s) = ^ /(27rm, 0.125a, (aV + (27r(m + l/2))2)V2) 

m=l 

+ 7(27rm, 0.125a, (aV + (27r(m- 1/2))^)^/^) (24) 
r°° X dx 

Jo (x2 + ng(x)/52)(x2 + + a2s2) ■ 



+ 



We now solve eqn (21) numerically. We do this by fixing a and N and varying 
s = Ti/a until the R.H.S is equal to unity. We note the value of s for which 
eqn (21) is satisfied for those values of and a. We then choose new values 
of a and N and repeat the process. If eqn (21) is not satisfied for any s at the 
values of a and chosen, the only solution to eqn (17) for constant mass is 
the trivial one s = 0. It is important to note that, in these calculations, the 
infinite sums are evahiatcd numerically using Mathematica: the convergence 
is very slow if an attempt is made to truncate them at a (large but) finite 
number of terms. 

In our analysis we also require the zero temperature limit of eqn (21) 
which is found to be 

i_ 1 r ^25) 

{2T:fN Jo (x2 + s2)(2; + 0.125) ^ ' 
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where x — \li\/a. On doing the integration on the R.H.S one obtains the 
following equation 



2 



V (0.125)2 + 5- 

inf^l. (26) 



(0.125)2 + s2 VO.125 

Now if s « 0.125 we get the approximate solution s ~ 0.25 exp(— 7r2A^/4), 
in agreement with [1], remembering that we are including only the longitudi- 
nal self-energy, which contributes one half of the total contribution at zero 
temperature, so that our "2A^" corresponds to Pisarski's 'W'.In fact, we 
have found that a more accurate fit to the exact numerical solution of (25) 
is provided by 

s ~ 0.25exp(-0.967rA^). (27) 



With the zero temperature result and our data for T 7^ 0, we are able to plot 
solutions of (21) as a function of 1/a = ksT/a for different values of N, as 
shown in fig. 4. It is important to see that unlike previous calculations our 
zero temperature results join on smoothly to our finite temperature results. 
This is because we have used a retarded form of the photon propagator 
with the correct T — > limit, and not the instantaneous form used in 
[9]. Another thing to note is that the shapes of the solutions as functions 
of ksT/a are markedly different from those of [9], [10] and [11]; they seem 
more to resemble the shape of the B.CS constant mass-gap solution. We 
also plot Tcks/cx against N, in fig.5. Tcks falls off exponentially and takes 
the following form, to very high accuracy: 

(rcA;B)/a ~ 0.081 exp(-0.967r7V). (28) 

^Prom eqn (27) and eqn (28) we can calculate r, where r — ^^^^^;we find 
r = 6.17. 



In this calculation there is therefore no Nc- The reason for this can be 
seen as follows. We first note that increasing s and decreasing a (increasing 
T) causes the R.H.S of eqn [21), -^Sl, to get smaller. The condition de- 
termining Nc is that \im.a^Qo{aSL) < 2tiNc for all s, since then eqn (21) can 
no longer be satisfied, and the only solution to the gap equation is s = 0. 
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Since 5^ is a monotonically decreasing function of s we need only show that 
Uma^oo(a.5'/,(s = 0)) is finite to estabUsh the existence of Nc- But one easily 
sees from eqn(26) that linia^oo (a5'L(s = 0)) is actually infinite, and so Nc 
does not exist. 

In [9] the instantaneous approximation was used with the constant gap 
approximation, giving a higher value r, r ~ 10. There is good agreement 
between our results and those of [12], in which retardation was included in a 
real time formalism. For example, in [12], using the "best average self-energy" 
(n^^ in the notation of [12]),s had the value 0.01205 for N — 2, correspond- 
ing to a value of s = 0.1225 from (27) with N = l;and {kBTc)/a was 0.0376, 
corresponding to 0.0397 from (28). These values lead to r = 6.41 in [12] as 
compared to our value of 6.17. The only other difference between our results 
and those of [12]-not a large one-is that we see very little change of r with 
A^,whereas [12] found a small variation with N . 

In the following section we shall look at solutions which arc frequency 
dependent and contain the retarded photon propagator calculated in section 
2. 

4 The fequency dependent gap equation 

In this section we depart from the constant gap approximation, and at- 
tempt to solve the retarded Schwinger- Dyson equation (17) allowing the gap 
function to depend on frequency only. Frequency dependence can be in- 
troduced relatively easily into the equation, unlike momentum dependence 
which leads to a highly non-trivial angular integral of the kernel of our equa- 
tion, which is unlikely to have an analytic result. The resulting computational 
complexity has deterred us from including momentum dependence.This is un- 
fortunate, for it would be interesting to compare results calculated with a 
full retarded 

momentum dependent S-D equation with the results of [10]. While this 
is true, introducing frequency dependence will give us new insight into the 
nature of the gap, and as in [10] we shall be able to calculate a value of N^. 
In the spirit of [10], 

we again neglect wavefunction renormalization. It has been pointed out 
in [3] that wavefunction renormalization should be introduced, but it has 
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been found in [11] that for a momentum dependent S-D gap equation, there 
still exists an Nc and r is httle affected. 



We start our study of the frequency dependent solution at eqns (17) and 
(18). We shall restrict ourselves again only to the longitudinal mode; this 
is because the i-r divergence in the tranverse mode is still present, although 
only in the mode for which go6 = ^o/ ~ Po/ = 0. We see that because E is 
now frequency dependent we can no longer remove a factor of S from each 
side of eqn (17), and A^^ is now dependent on po/ = ^(^ + 1/2) • But the 
angular integration is again trivial. On rearrangement we get: 

s{n)^^S-,{a,s) (29) 

where as before s = S/a and a = a(3. 5'£ is the frequency dependent 
longitudinal contribution, expressed as: 



S^ia, sin)) = I X dx 
Jo 



s{n) 



+ E 



[x^ + /?2n0)(a;2 + (27r(n + l/2))2 + a^s(n)) 
s{m) 

'^^ x^ + (27r(m - n))2 + 0.125a(a;2 + (27r(m - n))2)V2 

^ a;2+ (27r(m + l/2))2 + a2s2(m) 

where x = f3\h\.The first term represents the m = n mode, where qob = 0. 
Again it is possible to do the ^-integration. Using the functions, I{d, a, c) we 
are able to express 5'£ as: 

oo 

S1{a,s{n)) = 5^ s(m)7(27r|m-n|, 0.125a, (aV(m) + (7r(2m+l))^)^/^) 

r°° X dx s{n) 
^ Jo x^ + /32no x^ + (27r(n + l/2))2 + a2s2(n) ^ 

Eqn (29) represents an infinite system of non-linear coupled integral equa- 
tions, one for each mode, m. We can reduce the number of equations that 
need to be considered by the following particle-antiparticle symmetry prop- 
erty, 

s(m) = s{—m — 1) where m > 0. (32) 
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It is important that the R.H.S of eqn (32) is s{—m — 1) and not s{—m), 
since s(m) is a fermonic quantity; pof is related to m by po/ = ^("^ + 1/2)- 
Therefore m and —m — 1 represent pof and — po/ respectively. Using eqn (32) 
it is possible to cast 5'£"(a, s(n)) into a final form convenient for numerical 
calculations: 

oo 

Sl{a,s) = Yl s(m)[7(27r|m-n|, 0.125a, (a V(m) + ((7r(2m+l))2)i/2) 

+ (7(27r|m + n + 1|, 0.125a, (aV(m) + ((7r(2m + 1))^)^^^)] 

/■~ X dxsjn) 

Jo (x2 + /32n0(x))(x2 + (7r(2n+l))2 + a2s2(n) 



+ 

+ s(n)7(27r(2n + 1), 0.125a, (a^s^(n) + ((7r(2n + 1))^)^^^), 



the last term being the n = —m — 1 mode in eqn (31). Although eqn (33) 
looks cumbersome compared with eqn (31), the advantage is now that only 
positive values of n and m need be considered, thereby reducing the number 
of equations in the system (29). 

In our numerical calculation we will want to limit the number of equa- 
tions to a finite number, M. We can think of eqn (29) as a sort of non-hnear 
matrix problem, where instead of an infinite matrix we will restrict ourselves 
to an M X M matrix. The justification for doing this is simple: we expect our 
solution to fall off monotonically with n, and as n — > cx) we expect s{n) 0. 
The effective limit M on the matrix size is determined by the requirement 
that s{n — 1 = M) be one tenth the size of s(0), which is the maximum 
value of s as a function of n. In choosing this size of matrix, the contribution 
from modes for which n — 1 > M has a very small effect upon our solution. 
Note that whereas in the constant mass case (Section 3) we had to include 
all frequencies in our summations,when the mass is frequency- dependent a 
natural cut-off emerges. 

The method wc shall normally adopt ( "method I") in solving the equa- 
tions (29) will be to start with a trial function having the correct n ^ oo 
limit, namely s{n) — > 0. We insert this trial function into eqn(33). Then we 
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work out the resultant value of s(n) using eqn (29), and reinsert the result 
back into eqn (33). After several such iterations we check if our matrix size 
is sufficient, by seeing if s{n — 1) < O.ls(O); if it is appreciably larger we in- 
crease the size of our matrix. To do this, we increase the value of m at which 
we truncate the sum in eqn (33), and we increase the number of equations 
in eqn (29) we must solve. When we are satisfied with the size of our matrix 
we keep iterating until s(n) changes by less than 1% over ten iterations (or 
more depending on the rate of convergence to the exact solution). It should 
be appreciated that the non-linearity of the system of equations requires this 
further iteration, even after deciding on the size of the matrix. 

If the rate of convergence is very slow, which is the case near a phase 
transition, we use "method 11". In method 11 we choose two trial functions. 
One of them, si, is chosen to decrease in magnitude with each iteration us- 
ing method I; the other, S|, to increase in magnitude with each iteration. 
We then average the two functions together, s — 0.5{s^ + si). Then we use 
method 1 on s for ten iterations or more. If § increases in magnitude we set 
5|(n) = s{n); if s decreases in magnitude we set S|(n) = s(n), and we av- 
erage again. We repeat the process until S| and si vary little from each other. 

An important technical difficulty is that as T gets smaller the size of our 
matrix increases: as T — > 0,M — > cxo,and the larger M is, the slower the 
calculation of s{n). Fortunately most of the interesting behaviour lies in a 
region of temperature where the matrix size, M, is manageable. Although 
we cannot approach close to the T ^ limit from our finite temperature 
calculations, we can get around this difficulty by including an analytic T — > 
limit to eqn (28), 

«(^) = r dy s{y)my - ^, 1/8, ^y^ + dy)) 

a (zTrj^iV Jo a * 

+ {I{y + ^,l/8,^y^ + s(y))) (34) 

where y = kof/a. Again we have used the antiparticle-particle symmetry 
property of the gap, s(po/) = s(— po/), to reduce the number of equations 
we need consider. Although the range of integration is infinite, the kernel of 
eqn (34) provides an effective cutoff A, where A < 1. Therefore we need only 
consider s(^) from s(0) to s(A) in doing numerical calculations. A is chosen 
in the same way as M; namely at a point where s(A) is suffiently small so 
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that the contributions from pQf > A are neghgible. Our T 7^ results and 
our T — results are sufficient to provide all the information that we require. 

We show the results of both our T = and T 7^ calculations as func- 
tions of ^ for N=l in fig. 6. One can see that the T — curve joins smoothly 
on to the T 7^ solutions. One should observe that although we plot the 
T 7^ solutions as continuous curves, they arc actually discrete points with 
Pof = ^{n + 1/2). As functions of these functions fall off rapidly for 
^ < 0.1. This behaviour is roughly the same as was seen in momentum- 
dependent calculations [10], where for 0.02 < < 0.1 the solutions fell 
rapidly in — . It is this behaviour - the rapid fall in s( — ) for ^ < 0.1 
- which makes E(T — 0,po/ — 0) considerably smaller than E(T = 0) const, 
our result from the previous section, as shown in fig.7.E(T = 0,po/ = 0) is 
roughly the same as S(T = 0, p = 0) in [10]. This suggests that the constant 
mass -gap approximation is not a good one, except as an order of magnitude 
calculation for the quantity r. In fig. 8a we show the m = mode of these 
solutions for various N as functions of and in fig. 8b we show the m — 
mode of these solutions for various as functions of N. One notices that 

a 

the shape of the solutions in fig. 8a is the same as that of the constant mass- 
gap solution (fig. 4), which we have characterised as B.C.S -like. By looking 
at fig.8b one can sec that as well as a critical temperature, Tc, there is also a 
critical number, iVc. Nc is the value of A'^ above which s = for all T and po/, 
and this can be found from fig. 8b by locating the point where the gradient of 
the T = curve tends to infinity. We only need consider the m = mode at 
T = 0, since s(m) is a monotonically decreasing function of both frequency 
and temperature. So if s{T = 0,po/ = 0) = this condition must be true for 
any T and po/- We find that Nc hes in the range iV^ ~ 1.8 — 2. In fig. 9 we 
show a N-T phase diagram, where we have included the value Nc — 1.8. Here 
we see that the phase boundary agrees qualitatively with that calculated in 
[10]. 

The existence of Nc leads us to the conclusion that any frequency or mo- 
mentum dependence in the mass-gap equation induces a value of Nc- Our 
calculated value of Nc agrees well with [10], but [10] uses the instantaneous 
approximation and so has the incorrect T — limit. When compared with 
the correct T ^ limit of the 3-momentum -dependent mass-gap equation 
of [2], our value for Nc is smaller. (In [2] Nc was calculated to be 3.2). Prom 
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our results we are able to calculate a table of values for r, each at different 
values of A^, as shown in table. 1. These values are not much different from 
the value r ~ 6 calculated in section 2, but now there is a significant de- 
pendence of r on A^; the values of r are seen to fall with increasing N. The 
insensitivity of the r values to the introduction of a momentum dependence 
in the S-D equation was seen in [10], although for the instantaneous approx- 
imation; again this is roughly seen to be the case with the introduction of 
frequency dependence as shown by our results. 

This completes our analysis of the frequency dependent solutions. In the 
next section we discuss ways of handling the transverse mode and other ways 
of extending our calculation. 

5 Conclusion 

As we have seen already in section 4 the constant mass-gap approxi- 
mation is not a good one to make except possibly for the calculation of r, 
which is little changed by frequency dependence for the range of A'^-values 
considered here. Also we have seen that frequency dependence induces an 
Nc- Our result differs from [2], for we have been unable to introduce fre- 
quency and momentum dependence into our calculations while preserving a 
smooth T — i> limit, and also because we have no transverse contribution 
to the mass-gap equation. The reason why the introduction of momentum 
dependence is a hard problem at T 7^ is that three-dimensional Lorenz in- 
varance is lost, since frequency dependence becomes discretized (in the imag- 
inary time formalism). At T = one can do the angular integration with 
full 3-momentum dependence, because one can exploit the three-dimensional 
symmetry.At T 7^ this is not possible, due to the preferred (frequency) 
direction; instead we are faced with an integral in 2-dimensional momentum 
space, exactly the integral which is discussed at the begining of section 4. 

Up to now in our analysis we have chosen not to discuss the transverse 
contribution, and have neglected it from our calculations due to the i-r di- 
vergence in its zeroth mode. It is important to stress that the i-r divergence 
is present only when T 7^ 0. At T = 0, the effect of the transverse mode 
is merely to double the number of flavours, so giving us an Nc in the range 
3.8 — 4. To treat the i-r divergence at T 7^ one needs to regulate the 
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integral in eqn (20). A reasonable way to do this might be to introduce 
massive fermion propagators into our calculation of II^j,, the mass of the 
fermion propagators being calculated self consistently in eqn (19). A major 
problem with this refinement would be the requirement of self consistency 
in our calculation of the photon propagator when we have any frequency or 
momentum dependence in the fermion mass. Nevertheless, it might be worth 
considering the simpler problem of a constant mass in this sort of calculation; 
this may not give us reliable values of S(T = 0) and ksTc, but the value of r 
it gives might not change very much if frequency or momentum dependence 
were to be included. To tackle the harder problem of including frequency 
or momentum dependence in this extension to our work, we may be forced 
either to consider our calculation of 11^,^ as a separate integral equation, so 
that the calculations would involve two integral equations coupled togther, 
not one; or we may have to use an ansatz for our expected form of the fermion 
mass. 

As well as the above extension to our work involving the transverse contri- 
butions, another important step will be to introduce momentum dependence. 
Although we have stated that this is difficult due to the computational com- 
plexity, one may be able to approximate the kernel in such a way as to simplify 
the problem. One of the major goals in this series of numerical calculations 
for T 7^ should be to join on with the results of [2] in a smooth T — > 
limit, and this necessitates a successful treatment of momentum dependence. 

As to the question of QED^ as a model of superconductivity in the con- 
text of our calculations, consider the case N — l.If we crudely double the 
number of flavours so as to account for the transverse mode, this will effec- 
tively count as = 2 which is the value required for QED^ in a model of 
high Tc superconductivity [13]. We find that ksTc ~ 10~^q; in agreement 
with [10], once rescaled to agree with [2] as regards A'c.As pointed out in [10] 
this gives a value of a of the order of %eV for Tc ~ lOO^X.As stated in [10] 
this is still much higher than the typical Heisenberg exchange energies, but 
might be acceptable by suitable rescaling of the fermion lattice operators in 
the lattice model of [13]. It is interesting to note as a concluding remark 
that our suggested way to treat the i-r divergence in the zeroth mode of the 
transverse contribution may effectively reduce a by an order of magnitude 
or more,so on this basis alone it is a calculation worth considering. 
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Figure Captions 



Figure 1. 


Contributions to A^,^ to leading order in 1/A^ 


Figure 2 


(a). 


Hi as a function of explicit m and 3-momenta,p;,. 


Figure 2 


(b). 


112 as a function of explicit m and 3-momenta,p6. 


Figure 2 


(c). 


Ha as a function of explicit m and 3-momenta,pfe 


Figure 3 


(a). 


Top graph: The numerical data points and the approxima- 



tion f3piy/8 for the m ^ modes of Ili{m = l,pfe), shown for pi,f3 ranging from 
— 4. Bottom graph: The numerical data points and the approximation 
for the m 7^ modes of Ili{m = I,]?;,), shown for pi,P ranging from 
0-36. 

Figure 3 (b). Top graph: The numerical data points and the approxi- 
mation f3pb/8 for the m ^ modes of Il3{m = l,Pb), shown for pbf3 ranging 
from — 9. Bottom graph: The numerical data points and the approximation 
Ppb/8 for the m ^ modes of Uslm = l,Pb), shown for pbP ranging from 
0-36. 

Figure 4. The numerical results for the constant mass calculation as func- 
tions of for iV = 0.5, 1, 1.5. 

Figure 5. Plot of kBTja vs. N for the constant mass calculation; the 
points are the numerical values calculated from (21), the curve is an approx- 
imate form for Tc{N) as given by equation (28). 

Figure 6. The numerical results for the frequency dependent solution, 
E(po/, T) as a function of pof for various ksT/a a,t N — 1. 

Figure 7. A comparison of T^const, the mass for the constant mass calcu- 
lation ,with T^freq{m = 0,T = 0), the zero-frequency mass for the frequency 
dependent calculation, as functions of N. 

Figure 8 (a). The numerical results for the frequency dependent mass cal- 
culation as functions of 7^ for A'" = 0.5, 0.7, 1. 
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Figure 8 (b).The numerical results for the frequency dependent mass cal- 
culation as functions of for T/cs/o; = 0, 0.001, 0.002, 0.004, used to estimate 
Nc in the case of TkB/ot = 0. 

Figure 9. A phase diagram for the frequency dependent solution, showing 
the boundary between the massless and massive fermion phases in the N — T 
plane. 

Table Captions 

Table 1. The values of r calculated for various N in the frequency depen- 
dent mass calculation. 
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Fig. 9 
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